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Abstract. The global astrometric observations of a GAIA-like satellite were modeled within the PPN formulation 
of Post-Newtonian gravitation. An extensive experimental campaign based on realistic end-to-end simulations was 
conducted to establish the sensitivity of global astrometry to the PPN parameter 7, which measures the amount 
of space curvature produced by unit rest mass. The results show that, with just a few thousands of relatively 
bright, photometrically stable, and astrometrically well behaved single stars, among the ~ 10 9 objects that will 
be observed by GAIA, 7 can be estimated after 1 year of continuous observations with an accuracy of ~ 10 -5 
at the 3cr level. Extrapolation to the full 5-year mission of these results based on the scaling properties of the 
adjustment procedure utilized suggests that the accuracy of ~ 2 ■ 10 -7 , at the same 3cr level, can be reached with 
~ 10 6 single stars, again chosen as the most astrometrically stable among the millions available in the magnitude 
range V = 12 — 13. These accuracies compare quite favorably with recent findings of scalar-tensor cosmological 
models, which predict for 7 a present-time deviation, |1 — 7I, from the General Relativity value between 10 _J and 
10- 7 . 
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1. Introduction 

Since General Relativity's (GR) first appearance, 
many alternative gravity theories have been proposed. 
Experiments were then needed not only to test the validity 
of GR against Newton's theory, but also every alternative 
theory of gravity against all the others (Will ll99lfl l2001). 

The Parametrized Post-Newtonian formalism (PPN) 
takes the slow-motion, post-Newtonian limit of all the 
metric theories and exploits their similarities to give them 
a unique, coherent framework by introducing a set of 10 
PPN parameters. Each theory is then characterized by a 
particular value for each of the PPN parameters. 

The validity of the post-Newtonian expansion is con- 
fined to those physical situations in which the energy den- 
sity, regardless to its specific form, is small. Considering 
the various physical cases, this means that gravity fields 
(U), velocities («), pressure to matter density ratios (P/p), 
and internal specific energy densities II, are all small quan- 
tities so that (Ciufolini & Wheeler ggggl Will lTUMjl 
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where e <C 1 is the "smallness parameter". The slow- 
motion hypothesis also implies that, for any quantity A, 
its time derivative is much smaller than its space deriva- 
tive, so that (Misner et al.|1973) 

\dA/dt\ 

\dA/dx\ ~ e; 

this makes the PPN framework particularly convenient for 
Solar System experiments (where e.g. U/c 2 < 10 -5 ). Each 
experiment allows one to measure the values of some PPN 
parameters, so that the theories that do not match those 
values are ruled out. 

The most investigated among the PPN parameters are 
7 and f3, which measure respectively the amount of space- 
curvature produced by unit rest mass and the amount of 
non-linearity in the superposition law of gravity. 

In GR the parameters 7 and (3 are both equal to 1, 
while other theories aiming at the formulation of a uni- 
fied theory predict small deviations from the GR value. 
The most promising of such theories consider the exis- 
tence of a scalar field which, along with the usual metric 
tensor, participates in the generation of gravity. For this 
reason, such theories are usually called scalar-tensor the- 
ories. As a consequence, the value of 7 predicted by the 
scalar-tensor theories of gravity deviates from 1, and it 
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is generally an adjustable value. Moreover, Damour and 
Nordtvedt iTffSSl showed that scalar-tensor cosmological 
models contain a sort of attractor mechanism toward GR. 
In practice the value of 7 changes with time and, what- 
ever the value at the birth of the universe, it tends to 
1 as t — > 00. The present-time value of |1 — 7I (i.e. its 
deviation from GR) depends on the efficiency of the at- 
tractor mechanism, and Damour and Nordtvedt calculate 
a value between 10 -5 and 10 -7 . More recently, Damour et 
al. I|2002al l2002bjl have given a new estimation of |1 — 7I, 
within the framework compatible with string theory and 
modern cosmology, which basically confirms the previous 
results^ 1 ). 

Although optical astrometry provided the first direct 
test of GR by measuring the amount of light bending due 
to the gravitational pull of the Sun (Dyson, Eddington, 
and Davidson 1920), observing conditions from the ground 
have severely limited the possibility of this technique to 
contribute to the search for a scalar component of gravity. 
Radio astrometry has been able to achieve significantly 
better results, however it is again the observing condi- 
tions (especially the effects of ionosphere and troposphere) 
which are ultimately limiting its accuracy in measuring 
gravitational light deflection from the ground. Presently, 
7 is known to an accuracy of ~ 10~ 3 — 10~ 4 (Will 2 001(1 . 
therefore well above the level necessary to detect the scalar 
component predicted by Damour and Nordtvedt. 

The success of the mission Hipparcos has proven that 
going into space is the way for increasing the accuracy of 
astrometric measurements to more interesting values. By 
utilizing the measurements of Hipparcos, (Froeschle et al. 
1997 see also Mignard 2001J were able to derive 7 with 
an accuracy comparable to that of the best radio exper- 
iment (Lebach et al. 1995). Therefore, expectations are 
that the micro-arcsec (/^-arcsec) accuracy targeted by the 
space missions SIM, approved by NASA, and GAIA, ap- 
proved by ESA, will be able to measure 7 with a precision 
of 10- 5 -10" 6 (Danner and Unwin RMfy or better (ESA 
2000, de Felice et al.f2IQ2D|, positioning optical astrometry 
at the fore front of experimental gravitation. 

In this article we focus our attention on the GAIA 
project, a space astrometry mission that has recently been 
confirmed as a Cornerstone mission of the ESA program 
of scientific satellites. GAIA is expected to be launched 
not later than 2012, with a possible window of opportu- 
nity in 2010. For this reason ESA has devised a mission 
implementation plan that will be able to cope with a 2010 
launch. 

GAIA is a scanning astrometric satellite which builds 
on the successful concept of the Hipparcos mission, but 
with order-of-magnitude improvements in, e.g., the num- 



2«Ld, where a ha d is the 



1 Within this framework 1 — 7 
dilaton coupling to hadronic matter. Its value depends on the 
model taken for the inflation potential V(x) <x X™> X being the 
inflation. So the range of the expected deviations from GR is 
betw een ~ 3 .6- 10~ 5 (n = 4) and ~ 0.5- 10~ 7 (n = 2) (Damour 
et al. l2002b| . 



ber of objects observed and measurement precision. It 
will be essential for the scientific advancement in many 
branches of astronomy, and especially in the fields of stel- 
lar astrophysics and galactic astronomy (Perryman et al. 
2001). Here we report on the first thorough attempt at de- 
termining the accuracy with which GAIA could measure 
the parameter 7. This is done through a realistic end- 
to-end simulation which considers the satellite most rel- 
evant modes of observation, the relativistic environment 
in which such observations take place, and the expected 
single-measurement errors. 

2. The Parametrized Post-Newtonian (PPN) 
Model 

Using the Eddington-Robertson form of the PPN formal- 
ism (where only the parameters 7 e j3 are considered), we 
have developed a model based on the assumptions that 
the only source of gravity is a spherical and non-rotating 
Sun and the observer is moving on a circular orbit. This 
scenario is consistent with that discussed in our previ- 
ous papers (de Felice et al. IT53BI IM)T1 hereafter Paper I 
and Paper II, respectively), where a non-perturbative ap- 
proach was developed as a first attempt at modeling the 
GAIA observations in a rigorous relativistic environment. 
Therefore, the results discussed in this paper can directly 
be compared with the findings reported in those earlier 
works. 

Under the above hypotheses the PPN metric becomes 
(Misner et al., IW5jl 



ds 2 



M /M 
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dt 2 + 



l + 2 7 - 



M 



[dr 2 +r 2 (d9 2 + sin 2 d(j) 2 )} 



(1) 



where r, 9, and <fi are spherical coordinates centered on 
the Sun, t is coordinate time, and M is the geometrized 
mass of the Sun. 

As observable we consider the cosine of the angle -012 
between two stars; it can be expressed as 



cos ^12 = 



(2) 



where kf and k^ are the four-velocities of the photons 
reaching the observer (i.e. the four-vectors tangent to the 
photons' null-geodesies) and h a p — g a /3 + u a up is a tensor 
which projects to the rest frame of the observer, i.e., the 
GAIA satellite. Here g a @ and u a are the space-time metric 
and the observer's four-velocity, respectively. 

As expected, the expression we found for the four- 
velocities of the photons k a depends on both 7 and /3; 
then, they could be made part of the data reduction pro- 
cess as unknown parameters, in addition to those describ- 
ing position and velocity of the stars. However, (3 enters 
the k a as a second order term compared to 7 (see also 
Eq-JU); therefore, it was decided not to consider (3 in the 
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derivation of the observation equations (see below) but to 
set it to 1, namely to its value in GR. 

Following the method developed in Paper I and II, 
and after a lengthy calculation, we then obtained the lin- 
earized observation equation as function of the astromet- 
ric parameters (i.e., angular coordinates, parallaxes, and 
proper motions), and of 7 in the form ( Vecchiato l2UUfj|) 

-Bin Via (*) <%2 =A 1 8di(t )+B 1 6<p 1 (to) + 

Ci 5pi + Di 5iJ,e 1 + Ex + 

A 2 S9 2 (t ) + B 2 5<fo{t ) + 

C 2 5p 2 + D 2 5fie 2 + E 2 S^ 2 + 

FS 7 , (3) 

where the coefficients are derived from the differentiation 
of the right-hand side of Eq. (J2J with respect to each of the 
unknowns. 

As in Paper II, parallax and (angular) proper motions 
are defined as 

d(9 d<j> 
p^r e /r, „ = ^ = -; 

here r® is the Earth's mean orbital radius, r is the coordi- 
nate radial distance of the star from the Sun and t is the 
coordinate time. 

3. The end-to-end simulation 

The simulation follows the procedure used in the previous 
non-perturbative works (Paper I and II), the main change 
being the presence of the new unknown 7 that modifies 
the design matrix of the system of condition equations. 

First, we generate the set of true quantities, which de- 
fine initial location and temporal evolution of the stellar 
positions on the celestial sphere and the true value of 7. 
This was set to its GR value, i.e., 7=1. The corresponding 
catalog values are calculated, as usual, by adding suitable 
root-mean-square (rms) errors to the true values. Next, we 
find the stellar pairs which can be observed by a satellite 
that sweeps the sky following a Hipparcos-like scanning 
law. Once the stellar pairs are known the true angular dis- 
tances (from the true coordinates) are calculated, and the 
satellite observations are generated by perturbing these 
true arcs with the observational error in Tabled The cat- 
alog arcs (from the catalog coordinates) are also computed 
at this stage. Only arcs joining stars lying in different fields 
of view (FOVs) were counted without degeneration. 

The result of these two steps is the generation of the 
measured quantity, — smijj(t)5ip, and of the coefficients in 
Eq. (PJ for the construction of the linearized condition 
equations for all of the pairs observed during the mission 
lifetime we decided to simulate. Finally, the least-squares 
solution of the system is found by means of a conjugate- 
gradient method, suitable for large and sparse matrices 
like ours, and the errors are computed by direct compari- 
son to the true values. 



The least-squares solution returns the estimates of the 
adjustments to the catalog values of the unknowns in 
Eq. |(3J), i.e., SO, 5(f), 5[ig, dp, and £7. For the ad- 

justment to the 7 parameter, the corresponding rms error 
is indicated with the symbol as 7 - Notice that this error is 
the same as the error of 7, <r 7 , as 7 = j cat + 6~f, where j cat 
is the simulated catalog value of the deflection parameter. 
The symbols as-y and 07 are both used in the reminder of 
this article. 

4. Experiments and results 

4.1. An upper limit on the accuracy of "7 

Before going into the details of our experimental campaign 
it is useful to show, through a simple order-of-magnitude 
calculation, what kind of accuracy can be expected for 7 
with GAIA-like observations. This calculation starts from 
the consideration that the satellite measurements for the 
estimation of 7 can be thought of as "Eddington-like" 
measurements of very high precision but with the stars 
at some tens of degrees from the solar limb. Each ob- 
servation contributes to the determination of 7 with a 
precision of A 7 ~ 1CT 2 — 10~ 3 ( 2 \ If is the number 
of such observations, the final accuracy will be approxi- 
mately ct 7 ~ A 7 • N^ 1 / 2 . For a period of, say, 1.5 years of 
continuous operations and 6500 of the stars in the mag- 
nitude range V=12-13 mag, GAIA would provide (see be- 
low) about 450000 observations, thus cr 7 ~ 3 ■ 10~ 6 . In the 
calculation above we have disregarded the facts that the 
geometry of the GAIA observations is different (the grav- 
itational deflection is "seen" through its differential effect 
along the arcs joining the star pairs) and that 7 is not the 
only unknown to be estimated. Therefore, that value of 
<t 7 is clearly representative of a best case scenario, and it 
is used only for comparison in the following discussion. 

The simulation campaign was split in two parts. The 
first part was devoted to test the new code written for 
the PPN model and to make sure that the results were 
compatible with those in Paper I and II, and with the em- 
pirical prediction for 7 derived above. The second set of 
experiments was destined to establish the relation describ- 
ing the GAIA sensitivity to 7 as function of the number of 
stars. Such a relation can then be used to make predictions 
on the accuracy of 7 for any given number of stars. 

4.2. Validation of the PPN model 

Table ^ lists the values of the input parameters which 
were common to all of the simulation runs. The main dif- 
ference with our earlier work is in the single-observation 
error cr b s , which is set here to 10 /iarcsec. This is the value 

2 This can be seen by taking the expression for the light 
deflection of a light source (e.g. Misner et al.. 119731 . solving 
for 7 and then applying the error propagation formula. Then, 
a single observation of a star at 45° from the Sun with a b s = 
10 /xarcsec, yields A 7 ~ 2 • 10~ 3 . 
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Parameter 


Numerical value 


Comment 


orbital radius 


1.496 • 10 m 


same as Earth s orbital radius (it© J 


precession angle 


43° 


same as solar aspect angle 


precession speed of the spin axis 


6.4 rev/yr 




satellite spin period 


128 min 




angle between the 


54° 




viewing directions 






field-of-view 


1°6 




of each telescope 






mission starting time (to) 


-T/2 


minimum correlation between 






coords, and proper motions 


radius of the simulated sphere 


2 mas 


uniform density sphere 






of 500 pc in radius 


catalog error on 


2 mas 




coordinates and parallax 






catalog error on 7 


cr 7 = 2 ■ 10" 3 




single- measurement error (a bs) 


10 /larcsec 


as expected for pairs of 



V ~ 12 — 13 mag stars 



Table 1. Most relevant parameters common to all of the simulations. T is the mission duration. For the dynamical 
simulations, the stars are generated within a uniform density sphere of 500 pc centered on the Sun (i.e., p > 2 mas). 
Of course, the distribution of the resulting sample of simulated stars is representative only of a (relatively) small 
portion of the actual Galaxy. Nevertheless, as mentioned in Paper II, the adopted values are sufficiently realistic for 
the immediate scope of this work, which is to gauge the sensitivity of GAIA to the estimation of 7. 



expected for the error of one arc joining pairs of equal- 
magnitude stars, approximately 5 magnitudes brighter 
than those utilized in our previous experiments, i.e., V~12 
mag. The 10-/xarcsec error is compatible with current es- 
timates of the GAIA error budget for stars brighter than 

V = 12- 13 (ESA,nnoni) (3) - 

As for the 7 parameter, the starting (catalog) value for 
the PPN parameter is generated from the true value by 
adding an error of <t 7 = 2 • 10 -3 , which is comparable with 
current best estimates. 

The new code was tested on different sets of 50 simu- 
lations, and each Monte-Carlo set was run with the same 
values of the input parameters to have statistically signif- 
icant results. The results of this series of experiments is 
summarized in Table 

The first row is representative of the runs with a "static 
sphere" , i.e., the ideal situation where all the stars have no 
intrinsic motions and are located at such a large (infinite) 
distance that the parallactic motion due to the observer is 
also null. In this case, the location on the celestial sphere 
of the simulated stars is completely specified by their an- 
gular coordinates (6 and <fi). Therefore, each linearized ob- 
servation equation, Eq|3 has only five unknowns, the four 
corrections to the angular coordinates of each star pair, 
and the adjustment to the 7 parameter. 

The second row presents the results of the simulations 
with a "dynamical sphere", i.e., the case where the stars 

3 As mentioned in Paper I, it is assumed that the physical 
properties of the stars considered in these experiments are such 
that they do not show any intrinsic astrometric noise which 
adds to the measurement error. For example, single and non- 
variable solar-type dwarfs would be ideal targets. 



T 


n* 


^obs 


Q 


O~obs/ \fQ 


0~5p 






1.5 


6500 


453071 


34.85 


1.69 




1.76 


1.88 


2.0 


5000 


356061 


23.74 


2.05 


2.05 


1.31 


1.75 



Table 2. Summary of the results of the different sets of 
simulations utilized to test the new code developed for 
the PPN model. Q is the number-of-observations to the 
numbcr-of-unknowns ratio. 



move with time. For simplicity, only the parallax motion 
was simulated; the stars were generated within a uniform 
density sphere of 500 pc centered on the Sun, i.e., with 
parallaxes p > 2 mas (200 times the measurement error), 
and, as for the static case, with no intrinsic (cosmic) mo- 
tion. Notice that the mission duration was increased to 
2 years. According to the findings in Paper II, this is the 
minimum observation period required for a reliable recon- 
struction of the dynamical parameters. 

As expected, the true errors of the astrometric param- 
eters given in Table [2] compare quite well with what was 
obtained in Papers I and II after taking into account the 
differences in the values of Q and the factor of 10 in mea- 
surement errors, consistent with the 5-mag difference of 
the stars considered in the new experiments. In particu- 
lar, the results reported in the last row of Table 2 of Paper 
II, reproduce very closely what shown in our Table [5] for 
a mission duration of T=2 years. 

4.3. What accuracy on 7? 

Of much greater interest is the fact that the errors of 7 
after the sphere reconstructions shown in Table are con- 
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n* 


< nobs > 




5000 


1.78 ■ 10 b 


1.99 • 10 _b 


6000 


2.56 ■ 10 b 


1.91 • 10~ 5 


7000 


3.49 ■ 10 b 


1.23 ■ 10~ b 


8000 


4.56 ■ 10 s 


1.33 • 10 _b 


9000 


5.79 • 10 b 


1.34 • 10" 5 


10000 


7.14 ■ 10 b 


9.79 • 10 _b 


11000 


8.63 ■ 10 s 


9.19 ■ 10"" 


12000 


1.03 • 10 b 


1.07- 10~ b 


13000 


1.20 ■ 10" 


8.16 ■ 10" b 


14000 


1.40 • 10" 


7.87 ■ 10~ b 


15000 


1.60 ■ 10° 


6.19 ■ 10~ b 



Table 3. Computed errors on the estimation of the PPN 
parameter 7. The numbers refer to the results of the eleven 
T = 1 yr Monte-Carlo simulations with increasing n* 
(number of stars). In the second column < n bs > 1S the 
mean number of observations for the 50 simulations. 



GAIA-like astrometric mission . 

Number of observation r; ^ vs. number of si 




Fig. 1. Number of observations versus number of stars. 
The solid curve represents the interpolated law n b s =k ■ 
n* 2 , /c=0.0071264. 



siderably close to the best-case value of ~ 3- 10 6 derived 
at the onset of this section (see also Table 4 in de Felice 
et al 2000). 

This encouraging result brought us to consider a new 
set of simulations with the intent to study the accuracy 
of 7 for larger samples of stars. Indeed, the all-sky sur- 
vey nature of the GAIA observations ensures that all of 
the objects in the magnitude range of interest will be ob- 
served during the operational life of the satellite. This 
means that, although the actual number will considerably 
reduce because of the stringent requirements on the in- 
trinsic astrometric stability of the sources, millions of po- 
tential targets will be available for the "7 experiment". 
Unfortunately, the computing power needed to perform 
a data reduction simulation of such a size is beyond our 
present resources, therefore we had to resort to an alter- 
native schema to find the desired answer. We generated 
11 simulations of a static sphere, each consisting, as be- 
fore, of 50 runs with the same initial conditions and the 
mission duration set to 1 yr. The number of stars was in- 
creased from n* =5000 to the maximum extent possible, 
i.e., n*=15000. The results are listed in Tabled! 

A single simulation, of the 50 comprising each Monte- 
Carlo set, yields one value of the difference 67 = 7* — 
7(true 7 minus estimated 7); each simulation is then a 
measure of 7* affected only by random errors, and the 
sample of 50 S^s has a Gaussian distribution. The stan- 
dard deviation of this distribution, cr 7 , is the measure of 
the error on 7 of each Monte-Carlo group of simulations 
(Table 0J). For example, if a difference of 3 • 10~ 5 from 
the GR value, 7*=1, is measured for 7 with an error of 
07 ~ 6 • 10~ 6 , this would be interpreted as a 5-cr detection 
of a deviation from General Relativity. 

Fig-fflshows the relation n i, s — k ■ n* 2 we derived from 
interpolating the data in Table |21 Also, simple statistical 
considerations suggest to fit the results in columns 2 and 



3 of Table |21 to the ag-y oc " 
shown in Fig. [3 



-1/2 

obs 



relation, and the result is 




GAIA-Like astrometric mission 

results for the sensitivity on 5y 



MoiikvLLrlo results 
Fit wilhao. =kn for 



53073 
2e+07 



Fig. 2. crs-y versus number of observations. The solid curve 



represents the fit to the relation <t<5 7 
0.0088179 ±0.0000097. 



k • n 



-1/2 

'obs ' 



k 



We are now in a position to make realistic predictions 
on the error of 7 for a much larger number of stars (ob- 
servations), comparable to the size of the stellar sample 
expected to be surveyed by the satellite at the magnitude 
limit of interest here. For ~ 10 6 stars (~ 9.8 • 10 8 observa- 
tions), the value of ag-y estimated from the curve in Fig. [3 
is ~ 9 • 10~ 8 . This means that a 1 yr-long GAIA-like mis- 
sion could measure a value of 3 ■ 10~ 7 for |1 — 7], which 
would then represent a 3<7 detection of the deviation of 
the actual gravitation from GR. And at the end of the 
satellite lifetime, after 5 years of continuous observations, 
the sensitivity to a 3er detection will further improve by a 
factor of l/v2 to ~ 2 • 10~ 7 . These numbers imply that 
GAIA could reach the lower bound of the interval of possi- 
ble deviations predicted by Damour and Nordtvedt ( 1993) 
and by Damour et al. (2002b)l. 

4.4. How do we compare with the Hipparcos result? 

A Monte-Carlo experiment analogous to those above but 
with the relevant mission parameters set to the values of 
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the Hipparcos mission was also conducted. We deemed it 
as quite important to compare our findings to the work of 
Froeschle et al. (1997), who attempted the first direct de- 
termination of 7 by means of the global astrometric data 
taken with the Hipparcos satellite. Indeed, this would add 
confidence on the ability of our simulation to make realis- 
tic predictions on the possibility to derive a very accurate 
value for 7 with GAIA. 

As "reference" experiment we adopted the case with 
44000 stars in Table 1 of Froeschle et al. flMTj) . which 
resulted in an error on the deflection parameter of 07 = 
4- 10" 3 . 

We first simulated n*= 15000 stars and an observing 
period of 1 year; the observational and catalog errors were 
set to the values used in the Hipparcos experiments, i.e., 
Cobs = 3 milli — arcsec and cr ca t = 1 arcsec, respectively. 
The analysis of the usual 50 Monte-Carlo runs resulted in 
07 = 2 • 10 -3 . We then scaled this value to the duration of 
the Hipparcos mission, 3 yr, and to the number of obser- 
vations expected for 44000 stars (from the empirical law 
in Fig.l); this yielded 07 = 1 • 10~ 3 . 

That our experiment resulted in a much better value 
of 07 should not come as a surprise, simply because we 
put ourselves in a more favorable situation: the instru- 
ment and satellite attitude were both assumed perfect in 
the simulation. In particular, having assumed a perfect as- 
trometric instrument (optics, focal plane, and detectors) 
we disregarded any possible unmodeled systematic effect 
that could bias the estimation of 7, e.g., effects which 
would mimic a parallax zero-point error, as it has been 
the case for Hipparcos (Lindegren et al. 119921 Froeschle 
et al. I1997f) . In fact, in the observation equation derived 
from the reduction model used in the Hipparcos 7 experi- 
ment, one can see that the unknown which represents the 
parallax zero-point (common to all stars) is strongly cor- 
related with the 7 parameter. It can be shown that, for the 
Hipparcos mission parameters, such correlation amounts 
to p ~ —0.92 (Mignard 2001). A correlation of this mag- 
nitude, in turn, increases the error in the estimate of 7 by 
a factor of 1/ y/\ — p 2 ~ 2.6 ( - 4 - ) . Therefore, the simulated 
"replica" of the Hipparcos experiment is to be considered 
consistent within a factor of 1.5, and not a factor of 4, with 
the published results based on the real data; quite an en- 
couraging agreement given the quasi-ideal assumptions of 
our simulation. 



4 While it is extremely important that the astrometric pa- 
rameters be free of this kind of systematic effects which, though 
very small, could spoil any astrophysical result statistically in- 
ferred from such data, the problem of 7 is of completely differ- 
ent nature. In other words, we should be aware of systematic 
effects which correlate with the 7 parameter, but such effects 
need not necessarily be modeled, provided they are smaller 
than the deviation of 7 from the GR value we are trying to 
detect. 



5. Summary and conclusions 

The global astrometric observations of a GAIA-like satel- 
lite were modeled within the PPN formulation of Post- 
Newtonian gravitation. Although simplified (a spherical 
and non-rotating Sun is the only source of gravity), this 
PPN model has allowed, through extensive end-to-end 
simulations, a realistic evaluation of GAIA's sensitivity 
to the direct estimation of the light deflection parameter. 
The results show that the satellite could measure 7 to 
~ 10~ 7 (ler) after 5 years of continuous observations, and 
using a subset of approximately 10 6 stars chosen as the 
most astrometrically stable among the millions available 
in the magnitude range V = 12 — 13 of the GAIA survey. 
Notice that after just one full year of observations 7 could 
be estimated with an error only a factor of \[2 worse than 
the value above. 

A comparison with the Hipparcos results has provided 
a way to gauge the degradation factor to be expected in 
going from an ideal instrument to the real astrometric pay- 
load and satellite. The factor of 1.5 we found in that case 
is very encouraging; however, moving from the "milli-" to 
the "micro-" arcsec regime required by GAIA, the degrada- 
tion might become larger. Future work will have to deal 
with the fact that the accuracy of 10~ 7 sets a goal for 
both the observational model, which will have to include 
all the details and the implications of Solar System gravi- 
tation, and instrument development and modeling, which 
will have to concentrate in identifying and, possibly, re- 
move (through hardware improvements and/or calibration 
procedures), all relevant systematic error sources. 

This work has provided quantitative evidence that the 
micro-arcsecond global astrometry of GAIA appears ca- 
pable of testing general relativity to unprecedented levels. 
It will do so directly by accurately measuring the amount 
of light bending produced by gravity; a modern rendi- 
tion, about a century later, of the experiment of Dyson, 
Eddington, and Davidson, then the first proof of Einstein's 
theory. 
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